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A  coarse-grained  model  is  used  to  study  the  structure  and  dynamics  of  a  human  immunodeficiency 
virus- 1  protease  (1DIFA  dimer)  consisting  of  198  residues  in  an  effective  solvent  medium  on  a  cubic 
lattice  by  Monte  Carlo  simulations  for  a  range  of  interaction  strengths.  Energy  and  mobility  profiles 
of  residues  are  found  to  depend  on  the  interaction  strength  and  exhibit  remarkable  segmental 
symmetries  in  two  monomers.  Lowest  energy  residues  such  as  Arg41  and  Arg140  (most  electrostatic 
and  polar)  are  not  the  least  mobile;  despite  the  higher  energy,  the  hydrophobic  residues  (lie,  Leu, 
and  Val)  are  least  mobile  and  form  the  core  by  pinning  down  the  local  segments  for  the  globular 
structure.  Variations  in  the  gyration  radius  (. Rg )  and  energy  (Ec)  of  the  protein  show  nonmonotonic 
dependence  on  the  interaction  strength  with  the  smallest  Rg  around  the  largest  value  of  Ec.  Pinning 
of  the  conformations  by  the  hydrophobic  residues  at  high  interaction  strength  seems  to  provide  seed 
for  the  protein  chain  to  collapse.  ©  2010  American  Institute  of  Physics,  [doi:  10. 1063/1.3358340] 


The  catalytic  function  of  an  enzyme  such  as  1DIFA,  a 
human  immunodeficiency  virus  (HIV)-l  protease,  via  hy¬ 
drolysis  of  peptide  bonds1  is  believed  to  be  critical  for  the 
growth  of  HIV  in  host  cells.  Despite  a  considerable  interest 
in  recent  years,  the  search  for  an  ideal  inhibitor  is  far  from 
complete  due  to  mutation  of  HIV  protease.  Understanding 
the  structure  and  dynamics  of  such  a  protein  in  an  appropri¬ 
ate  solvent  is  therefore  very  important  in  the  search  for  pro¬ 
tease  inhibitors — a  targeted  effort  toward  finding  ways  to 
curtail  the  progression  of  HIV  infection.  Structural  stability 
of  HIV-1  protease  has  been  extensively  examined  using  ah 
atom6  molecular  dynamics  (MD)  and  coarse-grained9  MD 
simulations  in  an  effort  to  understand  the  effect  of  inhibitors. 
Most  of  the  analysis  seems  to  concentrate  around  an  astound- 
ingly  symmetric  structure  from  the  protein  data  bank  (PDB) 
involving  beta  sheets,  a  helices,  and  strands  resembling  the 
elements  of  a  face  involving  flaps  with  ears  and  eyes.1,6,9 
Recovering  such  a  structure  in  a  coarse-grained  model  is 
amazing  at  such  a  microscopic  detail,  i.e.,  a  specific  confor¬ 
mation  from  a  huge  ensemble  with  a  large  number  of  con¬ 
figurations  in  computer  simulation  models  where  invoking 
approximations  is  inevitable.  Of  course,  a  x-ray  snapshot  of 
the  structure  does  not  tell  much  about  ah  the  possible  con¬ 
formations,  nor  do  snapshots  from  computer  simulations 
guarantee  the  results  from  approximate  force  held,  con¬ 
straints,  or  ad  hoc  assumptions.  Yet,  comparison  of  the  quan¬ 
titative  results  of  computer  simulations9  with  the  experimen¬ 
tal  snaps  may  be  possible  but  somewhat  premature  or 
accidental.  It  is  not  feasible  to  explore  the  complete  confor¬ 
mational  phase  space  via  the  all-atom  MD  simulations.  All- 
atom  simulations  involve  some  level  of  coarse  graining  in 
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incorporating  the  effects  of  electrons  in  atoms.  Reproducing 
a  structure  in  any  simulation  (coarse-grained  or  all-atom)  is 
perhaps  accidental  or  a  result  of  constrained  exploration 
around  a  limited  phase  space.  We  constrain  our  goal  toward 
gaining  a  general  insight  and  predicting  qualitative  trends  in 
our  coarse-grained  model. 

The  last  few  decades  have  witnessed  an  enormous  surge 
in  interest10-29  in  modeling  proteins.  Much  attention  is  fo¬ 
cused  on  identifying  the  universal  characteristics,  e.g.,  fold¬ 
ing  pathways  via  analysis  of  the  energy  landscapes  of  protein 
chains  as  well  as  their  specific  characteristics  that  entail  local 
structures  to  understand  binding  to  pertinent  targets.  While 
all-atom  modeling  is  generally  appropriate  for  probing  local 
structures  in  short  time,  coarse  graining  is  almost  unavoid¬ 
able  to  address  large-scale  properties.  Using  a  coarse-grained 
Monte  Carlo  (MC)  simulation  we  have  recently  examined30 
both  local  and  global  properties  such  as  energy  and  mobility 
profiles,  structural  correlations  of  each  residue,  and  confor¬ 
mation  of  a  HIV-1  protease  (1DIFA  monomer)  which  con¬ 
sists  of  99  residues.  We  would  like  to  extend  this  study  to  a 
1DIFA  dimer  formed  by  combining  the  C-terminal  residues 
(96-99)  of  one  monomer  with  the  N-terminal  (residues  1-4) 
of  the  second  and  see  if  there  are  significant  changes  in  prop¬ 
erties. 


TABLE  I.  A  typical  interaction  matrix  among  H,  P,  E,  and  W  components. 
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I.  SEQUENCE  OF  1DIFA  DIMER  (HIV-1  PROTEASE) 

p1q2i3t4l5w6q7r8p9l10v11t12i13k14i15g16g17q18l19k20e21a22l23l24d25t26g27a28d29d30t31v32l33e34e35m36s37 

L38p39Q40j^41^Y42j^-43p44j^45jyj46j47Q48Q49j50Q5lQ52jr53j54j^-55y56j^57Q58Y59j)60Q61j62j^63j64£65j66^67Q68jj69j^70^71j72 
q^73  Y74y75j^76'y-77Q^78p79rp8Op81y82|^83j84j85Q^86p87|^88j^89j^9Orp9l092j93Q^94^95rp96j^97|^98p99plOOQ  101  jl02rp  103  j^104-yy  105 


QlO6plO7plO8LlO9yllOrpllljll2j^ll3jll4Qll5Qll6Qll7Lll8j^ll9£l2OYl2lLl22j^l23j)124rpl25Ql26Yl270l28j)129rpl3Oyl3lLl32 

£l33£l34jY|135gl36j^l37pl38Ql39pl40^yl41j^l42pl43j^l44jyjl45jl46Ql47Ql48jl49Ql50Ql52pl52jl53j^l54yl55pl56Ql57Yl58 

0l59Ql6Ojl6lLl62jl63£l64jl650l66Ql67jjl68j^-169Yl7Ojl7lQl72rpl73yl74pl75yl76Ql77pl78rpl79pl8Oyl81|^182jl83jl84Ql85 

pl86^187Ll88j^l89rpl9OQl91jl92Ql930l94rpl95Ll96p^l97pl98 


r 
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We  consider  the  coarse-grained  description  for  the 
1DIFA  dimer,  which  consists  of  198  residues  tethered  to- 
gether  on  a  cubic  lattice.  As  described  before,  the  empty 
lattice  sites  act  as  solvent  sites  for  effective  medium.  Each 
residue  interacts  via  a  Lennard-Jones  (LJ)  interaction 


uij=feij 


ra  <  r 


c 


with  surrounding  residues  and  solvent  (W)  sites  at  distance 
rtj  within  a  range  rc=  ^ 8;  cr=  1.  The  interaction  is  unique 
for  each  residue  type.  The  parameter/,  i.e.,  the  depth  of  the 
LJ  potential,  can  be  varied  to  control  the  quality  of  the  sol¬ 
vent.  The  characteristic  of  each  residue  is  captured  by  incor¬ 
porating  its  relative  hydrophobic  ( H ),  polar  ( P ),  and  electro¬ 
static  (E)  strengths  in  normalizing  their  interaction30 
presented  in  Table  I.  Interactions  (eEE)  among  the  electro¬ 
static  residues  Asp  {Ex),  Glu  (E2),  Lys  (E3),  and  Arg  ( E4 ) 

are!  £  1 1  =  £  1 2  =  £22  =  ^  33  =  ^  34  “  ^44  ~  ^  ■>  £13  —  £14— £23  —  £-24  ~ 

-0.4.  The  magnitude  of  each  interaction  element  (sij=£ji)  is 
based  on  the  insight  gained  from  the  all-atom  description  and 

30 

known  general  characteristics. 

Attempts  are  made  to  move  each  residue  randomly  with 
the  Metropolis  algorithm  subject  to  excluded  volume  con¬ 
straints  and  the  limits  on  changes  in  the  covalent  bond 
length.30  For  example,  a  randomly  selected  residue  at  a  site  i 
is  moved  to  one  of  its  randomly  selected  neighboring  lattice 
sites,  such  as  j  with  the  Boltzmann  probability  exp( 
-A Eij/T),  where  A £/■  is  the  change  in  energy  between  its 
new  (Ej)  and  old  (Et)  configuration  A Eij=Ej-Ei  and  T  is  the 
temperature  in  units  of  the  Boltzmann  constant  and  the  en- 
ergy  (fs  /y) .  The  attempt  to  move  each  residue  once  defines 
one  MC  step  time.  A  number  of  local  and  global  physical 
quantities  are  evaluated  during  the  simulations.  Physical 
quantities  in  our  study  include  energy  of  each  residue,  its 
mobility,  average  number  of  neighboring  residues,  and  mean 
square  displacement  of  the  center  of  mass  of  the  protein  and 
its  radius  of  gyration,  etc.  Simulations  are  performed  for  a 
sufficiently  long  time  (typically  for  1  X  107  time  steps)  with 
many  independent  samples  to  estimate  these  quantities  (see 
below).  Different  lattice  sizes  are  used  to  test  for  the  finite 
size  effects  on  the  qualitative  variations  in  these  physical 
quantities;  most  of  the  data  presented  here  are  generated  on  a 
3003  lattice. 


Figure  1  shows  the  snapshots  of  the  protease  at  the  end 
of  simulation  for  1  X  107  time  steps.  Small  and  large-scale 
globules  clearly  seen  in  this  figure  are  just  snapshots  among 
numerous  possible  conformations.  Visual  inspection  suggests 
that  the  hydrophobic  residues  constitute  the  core  of  these 
globule  formations.  The  conformational  phase  space  for  such 
a  protein  is  too  large  to  pin  down  a  unique  structure.  The 
global  quantities  such  as  the  radius  of  gyration,  energy,  trans¬ 
port  and  mobility  averaged  over  a  large  number  of  configu¬ 
rations  capture  their  general  characteristics.  Let  us  look  at 
some  of  the  local  quantities  such  as  energy  and  mobility 
profiles. 

The  energy  profile  of  each  residue  for  the  interaction 
strength  f- 150  and  200  is  presented  in  Fig.  2.  The  general 
pattern  in  residue  energy  appears  to  be  preserved  (with  / 
=  150,200)  with  somewhat  lower  energy  among  the  lowest 
energy  residues  at  higher  interaction  strength  (f=  200).  Least 
energetic  residues  are  easy  to  identify  which  includes 
Arg(R)41  and  Arg(/?)140,  the  residue  with  strongest  electro- 
static  strength.  Many  electrostatic  and  polar  residues  such 
as  Glu(E)21,  Glu(£)120,  Asp(D)25,  Asp(D)124,  Tyr( Y)59, 
Tyr(Y)158,  etc.  have  comparable  energy.  The  symmetry  in  the 
energy  of  the  residue  in  two  monomers  is  remarkable,  e.g., 
Pl-E of  the  first  monomer  and  P100-F198  of  the  second  such 
as  residues  R 41  and  Ruo ,  F59  and  F158,  etc.  Perhaps  this  is  the 
cause  for  the  symmetric  structure  for  the  protease  with 
flaps1,6,9  pointed  out  in  the  beginning.  Based  on  the  energy 
profiles  one  can  identify  the  segmental  stability  or  instability. 
For  example,  it  is  easier  to  find  the  least  energy  segment 
(D29-E35  ,Dm-E134),  the  most  energetic  segments 

\ 

FIG.  1.  [pink  (hydrophobic,  H),  gold  (polar,  P),  blue  (electrostatic,  £)] 
Snapshots  of  the  1DIFA  dimer  at  the  end  of  1  X  107  time  steps  with  the 
interaction  strength /=  150,  200,  and  250  on  a  3003  lattice. 
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FIG.  2.  Energy  of  each  residue  in  HIV-1  protease  (1DIFA  dimer)  with  198  residues;  two  monomers  (1-99  and  100-198  in  the  sequence)  are  separated  for 
clarity.  Sample  size  3003  is  used  for  time  steps  107  with  50  independent  runs  for  each  interaction /=  150  and  200. 


(A71-V82,A170-V181),  and  so  on.  Thus,  the  energy  profiles  pro¬ 
vide  a  useful  map  of  the  segmental  energy  and  possibly  their 
reactivity  with  regard  to  proteolysis. 

The  corresponding  mobility  profile  of  the  residues  is  pre¬ 
sented  in  Fig.  3.  The  mobility  (Mn)  of  each  residue  is  defined 
as  the  average  distance  traveled  per  unit  time  step.  Obviously 
the  residues  are  more  mobile  with  lower  interaction  strength 
(f=  150)  than  those  with  higher  interaction  strength  (f 
=  200). 


The  symmetry  in  the  mobility  profiles  of  the  residues  in 
two  monomers  in  1DIFA  dimer  persists  as  with  their  energy 
(Fig.  2).  The  distinction  among  different  segments  of  the 
protein  based  on  residue  mobility  becomes  more  transparent 
at  higher  interaction  strength  (f=  200).  The  most  mobile 
longest  segment  Asp25-Lys45  in  the  first  monomer  and 
Asp124- Asp 144  in  the  second  monomer  consists  of  polar  and 
electrostatic  residues  with  relatively  smaller  numbers  of  hy¬ 
drophobic  residues  (Gly,  Ala,  and  Met).  Some  of  these  highly 


FIG.  3.  Mobility  of  each  residue  in  HIV-1  protease  (1DIFA  dimer)  with  198  residues;  two  monomers  are  separated  for  clarity.  Sample  size  3003  is  used  for 
time  steps  107  with  50  independent  runs  for  each  interaction /=  150  and  200. 
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FIG.  4.  Root  mean  square  displacement  of  the  center  of  mass  of  the  protein 
(1DIFA  dimer)  with  the  time  step  for  various  interaction  strength  parameters 
/=  150-260;  the  slopes  of  the  long  time  data  sets  for  two  interaction 
strengths  (f=  150,260)  are  indicated.  The  inset  shows  the  evolution  of  the 
radius  of  gyration  (R  )  with  the  time  steps  for  representative  values  of  / 
=  150,  200,  and  260.  Larger  numbers  (1000-2000)  of  independent  samples 
were  used  for  data  sets  with  shorter  time  (up  to  105)  steps  and  lower  num¬ 
bers  (50-80)  of  samples  for  long  time  data  sets  on  a  3003  lattice. 


mobile  residues  are  also  the  least  energetic 
(Asp29-Glu35,  Asp128-Glu134)  which  are  interspersed  between 
high  mobile  segments  Thr26-Ala28  in  the  first  monomer  and 
Thr125-Ala127  in  the  second  and  Ser37-Pro44  in  the  first  mono¬ 
mer  and  Ser136-Pro143  in  the  second.  There  are  several  highly 
mobile  small  segments  separated  by  least  mobile  residues 
and  segments.  For  example,  in  the  first  monomer,  these  high 
mobility  segments  are  Trp6-Pro9,  Thr12,  Gly16-Gln18, 
Arg57-Gln61,  and  Gly78-Pro81  which  are  separated  by  least 
mobile  residues  and  segments  lie3,  Leu5,  Leu10,  Val11,  lie13, 
lie15,  Ile62-Ile64,  lie66,  Cys67,  lie72,  Val75-Val77,  Val82,  He84, 
lie85,  Leu89,  Leu90,  lie93,  Cys95,  and  Leu97.  The  same  residues 
are  also  in  the  second  monomer.  Note  that  all  the  least  mo¬ 
bile  residues  are  hydrophobic  which  seem  to  pin  down  the 
movements  of  neighboring  residues.  The  mobility  profile  de¬ 
pends  on  the  sequence  and  provides  a  road  map  of  specificity 
of  the  enzyme.  Despite  symmetry  in  energy  and  mobility 
profiles  in  two  monomers  of  the  dimer,  there  are  some  asym¬ 
metries  (e.g.,  Pl).  Regardless  of  the  remarkable  symmetries, 
it  is  difficult  to  identify  the  configurations  resembling  a  face 
with  flaps.1,6,9 

How  does  the  mobility  of  individual  residues  collec¬ 
tively  affect  the  global  dynamics  of  the  protein?  Variation  in 
the  rms  displacement  (. Rc )  of  the  center  of  mass  of  the  pro¬ 
tein  with  the  time  steps  ( t )  as  presented  in  Fig.  4  for  a  range 
of  interaction  strength  (f=  150,260)  can  provide  some  in¬ 
sight  into  the  global  dynamics.  The  slope  of  Rc  versus  t  on  a 
log-log  scale,  i.e.,  the  power-law  exponent  v  in  Rc^tv,  de¬ 
scribes  the  characteristic  dynamics.  At  a  relatively  low  inter¬ 
action  strength  (f=  150),  the  increase  in  Rc  with  the  time  step 
with  a  relatively  high  rate  shows  that  the  protein  chain  con¬ 
tinues  to  diffuse  (v~  1/2)  freely  in  the  long  time  regime. 
This  is  consistent  with  the  relatively  high  mobility  profile  of 
the  residues  (see  Fig.  3  with  /=  150).  The  dynamics  of  the 
protein  slow  down  on  increasing  the  interaction  strength  es¬ 
pecially  in  the  asymptotic  regime  after  its  relaxation.  With 
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FIG.  5.  Radius  of  gyration  (Rg)  left  and  energy  (Ec)  of  1DIFA  dimer  vs  the 
interaction  strength  /  in  steady  state  with  50  independent  samples  on  a  3003 
lattice. 


the  strong  interaction  strength  (/*>  220),  the  global  dynamics 
become  very  slow  (z7~0.14)  as  the  hydrophobic  residues  pin 
down  the  protein  chain;  the  onset  of  pinning  is  already  seen 
in  the  mobility  profile  of  the  residues  with  /=  200  (Fig.  3). 
The  quality  of  solvent,  a  measure  of  the  interaction  strength, 
thus  dictates  the  global  dynamics  of  the  protein. 

How  does  the  global  conformation  of  the  protein  relax 
starting  from  an  initial  random  configuration?  The  inset  in 
Fig.  4  shows  the  evolution  of  the  radius  of  gyration  (. Rg )  of 
the  protein  with  the  time  steps  for  representative  values  of 
the  interaction  strength,  i.e.,  /=  150,  200,  and  260.  Note  that 
a  million  time  steps  are  not  long  enough  to  reach  a  steady- 
state  value  with  f- 150,200  and  therefore  simulations  had  to 
be  continued  for  much  longer  time  (1  X  107)  steps.  The 
asymptotic  value  of  Rg  exhibits  large  fluctuations  in  the  long 
time  regime  (f~  106-107  steps);  the  statistical  fluctuations 
can  be  reduced  by  increasing  the  number  of  independent 
samples  (compare  short  and  long  time  data  sets  in  the  inset 
Fig.  4).  The  steady-state  value  of  Rg  is  evaluated  by  averag¬ 
ing  data  in  the  asymptotic  regime,  i.e.,  f~  106-107  for  Rg 
with  f=  150  and  5  X  106-107  for  Rg  with  /=  200.  Strong 
interactions  (f=  260)  pin  down  the  conformations  to  their 
smaller  values  of  Rg. 

Variations  in  the  steady- state  gyration  radius  (. Rg )  and 
energy  (Ec)  of  the  protein  with  the  interaction  strength  (f)  are 
presented  in  Fig.  5.  Despite  large  error  bars,  both  Rg  and  Ec 
show  nonmonotonic  dependence  on  the  interaction  strength 
with  opposite  trends.  The  smallest  size  of  the  protein  (. Rg  in 
left  figure)  is  not  near  the  minimum  energy  (right  figure)  of 
the  protein  around  f- 130.  In  fact,  the  smallest  value  of  the 
gyration  radius  appears  around  the  interaction  strength  / 
=  230  where  the  protein  energy  is  closer  to  its  maximum 
value.  Pinning  of  the  conformations  by  the  hydrophobic  resi¬ 
dues  at  such  high  interaction  strengths  seems  to  provide  seed 
for  the  protein  chain  to  collapse  (see  Figs.  1  and  3).  Driving 
the  protein  toward  its  globular  structure  by  increasing  the 
interaction  strength  leads  to  higher  energy  of  the  protein,  a 
somewhat  contradictory  observation  to  the  general  notion  of 
funneling  pathways  in  the  complex  energy  landscape.  These 
observations  for  both  local  and  global  physical  quantities  for 
the  dimer  protein  are  similar  to  what  we  have  observed  for 
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the  monomer.  Two  monomers  in  the  1DIFA  dimer  simply 
help  enhance  the  active  sites  (energy,  mobility,  and  pinning). 

Regardless  of  the  remarkable  symmetries  in  these  char¬ 
acteristics  of  the  dimer,  we  are  unable  to  see  the  signal  of  the 
stable  structure  that  resembles  a  face  with  flaps.1,6,9  It  is  dif¬ 
ficult  to  search  all  conformational  phase  spaces  in  such  a 
large  ensemble  of  the  protein  structure.  The  importance  of 
sampling  used  in  our  MC  simulation  via  the  Metropolis  al¬ 
gorithm  provides  an  efficient  searching  procedure  for  the 
most  accessible  states  (conformations)  energetically.  How¬ 
ever,  it  does  not  capture  all  states.  Therefore,  it  is  not  feasible 
to  look  at  each  configuration  and  identify  the  precise  folded 
configuration  explicitly  as  provided  in  the  PDB.  Apart  from 
coarse  graining  the  atomistic  details  the  selective  sampling 
though  efficient  and  feasible  may  be  one  of  the  unavoidable 
limitations  of  the  MC  simulations.  Such  an  approach  is  more 
valuable  in  identifying  the  trends  in  the  ensemble  averaged 
properties,  e.g.,  correlation  between  the  rapid  decay  of  the 
radius  of  gyration  with  the  interaction  and  the  energy  (Fig.  5) 
than  identifying  the  specific  structure.  One  of  the  primary 
implications  of  our  study  is  that  one  should  not  rely  on  a 
single  configuration  to  interpret  experimental  observables. 
Instead,  one  should  consider  ensemble  averaged  quantities 
for  trends  in  observables  including  relaxation  from  segmen¬ 
tal  characteristics  to  global  structure. 

In  conclusion,  MC  simulations  are  performed  to  examine 
the  local  and  global  properties  of  a  HIV-1  protease,  1DIFA 
dimer  consisting  of  198  residues.  A  coarse-grained  descrip¬ 
tion  is  used  to  model  the  protease  on  a  cubic  lattice  where 
empty  lattice  sites  constitute  the  effective  solvent  medium.  A 
residue  is  represented  by  a  node  of  the  protease  chain  and 
interacts  with  other  residues  within  the  chain  and  with  the 
solvent  sites  within  a  range  of  interactions.  The  interaction 
strength  is  varied  via  a  parameter  /  that  represents  the  depth 
of  the  LJ  potential  well.  Both  local  and  global  properties, 
conformation  and  dynamics  depend  on  the  interaction 
strength.  The  energy  profile  of  the  residues  in  the  steady  state 
reveals  that  the  segments  of  proteins  that  consist  of  polar  and 
electrostatic  residues  equilibrate  to  their  lowest  energy  con¬ 
figuration.  For  example,  some  of  the  low  energy  residues  are 
Glu21,  Asp  CD)25,  Arg41,  and  Tyr(F)59  in  first  monomer  and 
corresponding  residues  in  the  second  monomer  are 
Glu(£)120,  Asp(D)124,  Arg(R)140,  and  Tyr(F)158.  Hydropho¬ 
bic  residues  (lie,  Leu,  and  Val)  possess  higher  energy  but  are 
less  mobile  than  their  neighboring  electrostatic  and  polar 
groups.  Local  segmental  dynamics  and  structure  are  pinned 
down  by  these  hydrophobic  residues,  which  constitute  the 
core  of  the  globules  at  higher  interaction  strength.  There  is  a 
remarkable  symmetry  in  both  energy  and  mobility  profiles  of 
the  residues  in  each  monomer  of  the  1DIFA  dimer.  Despite 
these  symmetries  in  energy  and  mobility  profiles  of  each 
monomer,  it  is  difficult  to  identify  the  symmetric  configura¬ 
tions  resembling  a  face  with  flaps  due  to  large  phase  space 
with  enormous  numbers  of  configurations. 

Variation  in  the  RMS  displacement  with  the  time  steps 


shows  that  the  protein  diffuses  in  a  solvent  with  relatively 
weaker  interaction.  The  global  motion  of  the  protein  slows 
down  systematically  on  increasing  the  interaction  strength, 
where  pinning  by  the  hydrophobic  core  residues  induce  glob¬ 
ule  formation.  Variations  in  the  gyration  radius  ( Rg )  and  en¬ 
ergy  (Ec)  of  the  protein  show  non-monotonic  dependence  on 
the  interaction  strength  with  opposite  trends.  The  smallest 
value  of  the  gyration  radius  appears  around  the  interaction 
strength  (f=  230)  where  the  protein  energy  is  toward  its 
maximum  value.  Pinning  of  the  conformations  by  the  hydro- 
phobic  residues  at  such  high  interaction  strength  seems  to 
provide  seed  for  the  protein  chain  to  collapse. 
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